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We study the electronic states for Ca2- a; Sr a ;Ru04 in 0.5 < x < 2 within the Gutzwiller approx- 
imation (GA) on the basis of the three-orbital Hubbard model for the Ru t2 S orbitals. The main 
effects of the Ca-substitution are taken into account as the changes of the dp hybridizations between 
the Ru 4d and O 2p orbitals. Using the numerical minimization of the energy obtained in the 
GA, we obtain the renormalization factor (RF) of the kinetic energy and total RF, which estimates 
the inverse of the mass enhancement, for three cases with the effective models of x — 2 and 0.5 
and a special model. We find that the inverse of the total RF becomes the largest for the case of 
x — 0.5, and that the van Hove singularity, which is located on (below) the Fermi level for the special 
model (the effective model of x = 0.5), plays a secondary role in enhancing the effective mass. Our 
calculation suggests that the heavy fermion behavior around x = 0.5 comes from the cooperative 
effects between moderately strong Coulomb interaction compared to the total bandwidth and the 
modification of the electronic structures due to the rotation of RuOe octahedra (i.e., the variation 
of the dpn hybridizations and the downward shift for the d xy orbital). We propose that moderately 
strong electron correlation and the orbital-dependent modifications of the electronic structures due 
to the lattice distortions play important roles in the electronic states for Ca2- a; Sr a ;Ru04. 

PACS numbers: 71.27.+a, 74.70.Pq, 



I. INTRODUCTION 

Strongly correlated electron systems with orbital de- 
grees of freedom have attracted much interests due to a 
diversity of phenomena, e.g., a structure-sensitive super- 
conducting transition temperature,— a colossal magneto- 
resistance)^ and an orbital ordering and the correspond- 
ing collective excitations Among them, the layered 
perovskite Ca2- a; Sr a ;Ru04 has rich ground states^i 7 - in 
which the Ru Ad orbitals play important roles in deter- 
mining the electronic states. 

In Ca2- a; Sr a ;Ru04, substitution of Ca 2+ for Sr 2+ 
causes three kinds of lattice distortions, i.e., rotation, 
tilting, and flattening, which affect the electronic struc- 
tures^ Starting from x = 2 (i.e., Sr 2 RuO,d), which is a 
spin-triplet superconductor^— the ground state changes 
from the spin-triplet superconductor to a paramagnetic 
(PM) metal in the range of 0.5 < x < 2. In partic- 
ular, in the range of 0.5 < x < 1.5, the Ca substitu- 
tion induces a rotation of RuOg octahedra around the c 
axis (see Fig. [I}, while the crystalline structure remains 
tetragonal with a unit cell doubled. Except the vicinity 
of x = 2, the spin susceptibility behaves like the Curie- 
Weiss type with a Curie constant corresponding to nearly 
S = 1/2, and the Weiss temperature reaches nearly zero 
at x = 0.5£ In the range of 0.2 < x < 0.5, there oc- 
curs a structural transition from tetragonal phase to or- 
thorhombic phase at low temperatures, and the Ca sub- 
stitution induces the tilting of RuOg octahedra around 
a non-symmetric axis in the basal plane. Also, in this 
range, the spin susceptibility has a peak at some temper- 
ature, and the ground state becomes an antiferromagnet- 
ically correlated metali 7 - In the range of < x < 0.2, a 
flattening of RuOg octahedra along the c axis is induced 



and there is a first-order metal-insulator transition: the 
ground state becomes the antiferromagnetic insulator^ 

Most surprisingly, the coefficient of the electronic spe- 
cific heat, 7o, at low temperatures becomes large around 
x = 0.5. The largest value of j c reaches 255 mJ/mol- 
Ru K 2 at x = 0.5; this value is about 25 times of 
that obtained in the density-functional calculation for 
Sr2Ru04 within local-density approximation (LDA).— A 
similar mass enhancement has been observed in the op- 
tical conductivity measurement by using the extended 
Drude modelJ^ In addition, the spin susceptibility and 
Hall coefficient measured by applying the magnetic field 
perpendicular to ab plane show peak structures at T* ~ 
10 K for x = 0.3,i& which arc similar to those obtained for 
UPt3ii£ These experimental results indicate a formation 
of heavy fermions (HFs) around x = 0.5. 

Although there have been extensive theoretical works, 
the origin of this HF behavior has not been clarified yet. 
Some studie s 17 ' 18 based on the dynamical mean-field the- 
ory (DMFT) have proposed that the HF behavior results 
from the orbital-selective Mott transition (OSMT) for the 
d xz and dy Z orbitals; the OSMT is defined as a partial 
disappearance of quasiparticles only for some of the con- 
ducting bands. However, this contradicts with another 
DMFT calculation-^ which has used a realistic model of 
Ca 2 -xSr x Ru0 4 in 0.5 < x < 2: the OSMT does not 
appear. 

The OSMT for the d xz and d yz orbitals also contra- 
dicts with several experiments in Ca2-xSr K Ru04 around 
x = 0.5. The angle-resolved photoemission spectroscopy 
(ARPES) measurement 2 ^ for x = 0.5 and 2 has shown 
that the Ca substitution does not modify the topologies 
of the Fermi surfaces (FSs) for the d xz and d yz orbitals. 
This result is inconsistent with the proposal of the OSMT 



since the partial Mott gap should drastically affect the 
topologies of the FSs. In addition, the optical conductiv- 
ity measurement*^ for Ca2- a; Sr a ;Ru04 in 0.06 < x < 2 
has found that the Drude weights depend weakly on the 
Ca concentration. This is consistent with the ARPES 
measurement since the latter shows a small change of 
the FS from that for Sr 2 Ru0 4 . Note that the FS for the 
d xy orbital changes from an electron pocket for x = 2 to 
a hole pocket for x = 0.5. 

The aim of this paper is to discuss the qualitative 
origin of the HF behavior in particular near x = 0.5 
for Ca2- K Sr a ;Ru04. Wc study the electronic states for 
0.5 < x < 2 within the Gutzwiller approximation (GA) 
on the basis of the three-orbital Hubbard model for the 
Ru t2 g orbitals (i.e., d xz , d yz , and d xy orbitals). Wc 
assume that the Ca substitution affects the electronic 
structures mainly by the changes of the dp hybridiza- 
tions between the Ru 4<i and O 2p orbitals. Actually, the 
density-functional calculation within the LDA has found 
that the rotation of RuOg octahedra affects the electronic 
structures for Ca 2 - x Sr ;E Ru04 in 0.5 < x < 2.— The GA 
is used to include the effects of electron correlation non- 
perturbatively, in which the effects lead to the renormal- 
ization of the kinetic energy32r— 

Using the numerical minimization of the energy ob- 
tained in the GA, we obtain the renormalization factor 
(RF) of the kinetic energy for the Ru t^g orbtials and 
total RF, which estimates the inverse of the mass en- 
hancement, for the effective models of x = 2 and 0.5. 
We find that the difference between the total RFs for the 
cases of x = 2 and 0.5 becomes large as the intraorbital 
Coulomb interaction is strong; this arises from the criti- 
cality approaching the usual Mott transition, where the 
occupation numbers for the d xz / yz and d xy orbitals are 
1 and 2, respectively. In addition, we analyze a special 
model in which the van Hove singularity (vHs) for the d xy 
orbital is located on the Fermi level. We find that the 
total RF is smaller for the effective model of x = 0.5 than 
that for the special model; the vHs plays a secondary role 
in enhancing the effective mass. These results are con- 
sistent with the experimentally observed tendency of the 
effective mass in 0.5 < x < 2, i.e., monotonic increase of 
7e towards x = 0.5 Our calculation suggests that the 
HF behavior around x = 0.5 comes from the cooperative 
effects between moderately strong Coulomb interaction 
compared with the total bandwidth and the significant 
modification of the electronic structures for the Ru tig 
orbitals due to the rotation of Ru06 octahedra; the lat- 
ter includes both the variation of the dpir hybridizations 
and the downward shift for the d xy orbital. 

The paper is organized as follows. Section II is de- 
voted to the explanations of the method to take account 
of the main effects of the Ca substitution on the elec- 
tronic structures in 0.5 < x < 2 and the GA for a PM 
state of the degenerate d xz and d yz orbitals and the d xy 
orbital. In Sec. Ill, we show the numerical results of the 
GA for three cases with the effective models of x = 2 and 
0.5 and the special model. In Sec. IV, we compare our 
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FIG. 1: Schematic pictures of a xy plane of RuOe octahedra 
for (a) 1.5 < x < 2 and (b) 0.5 < x < 1.5. Black (white) 
circles represent Ru (O) ions. (f> is the angle of the rotation of 
RuC>6 octahedra. A and B are the indices of two sublattices. 
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FIG. 2: (Color online) Band structure of the Ru t2 g orbitals 
for the effective model of x = 2. The dashed black line repre- 
sents the chemical potential. 

results with previous theoretical studies and remark on 
the correspondence of our results with experimental re- 
sults. The paper concludes with a summary of our results 
in Sec. V. 



II. FORMULATION 

In the following, we choose the coordinates, x, y, and 
z, in the directions of the bonds between Ru and O ions 
at tfi — 0°, where <fi is an angle of the rotation of RuOe 
octahedra [see Fig. HR)]. Namely, the coordinates are 
fixed in the directions of Sr2Ru04 with = 0° even for 
the cases with finite values of 4>. For convenience, the 
five Ru 4c? orbitals, d xz , d yz , d xy , d x 2_ y 2, and d 3z 2_ r 2, are 
labeled 1, 2, 3, 4, and 5, respectively. 

According to the density-functional calculation o 27 ' 28 
for Sr2Ru04 within the LDA, the antibonding bands of 
the Ru t 2g and O 2p orbitals form the conducting bands 
in the vicinity of the Fermi level, and the density of states 
(DOS) near the Fermi level is originated mainly from the 
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FIG. 3: FSs for the effective model of x = 2. 
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FIG. 4: (Color online) DOS for the effective model of x = 2. 
The dashed black line represents the chemical potential. 



Ru t2 g orbitals. The partial density of states (pDOS) for 
the O 2p orbitals is roughly a quarter of those for the 
Ru tig orbitals. In addition, the x-ray measurements in 
Ca2_ 2; Sr :Z ;Ru04 have shown that the crystalline-electric- 
ficld (CEF) energy between the Ru e g orbitals and Ru t^g 
orbitals is of the order of 1 eVi 29 ' 30 Thus, the Ru t-i g or- 
bitals play main roles in determining the electronic states 
for Ca2- a; Sr a ;Ru04 at low temperatures. In Sec. IV, we 
remark on the roles of the Ru e g and O 2p orbitals. 

In order to study the electronic states, we use the three- 
orbital Hubbard model for the Ru ti g orbitals. The unit 
cell contains one Ru atom for 1.5 < x < 2, while it con- 
tains two Ru atoms for 0.5 < x < 1.5 due to the rotation 
of RuC>6 octahedra. As shown in Fig. [T] (b), the di- 
rections of the rotation alternate in the two-dimensional 
square lattice, leading to a unit cell doubled. 

The noninteracting Hamiltonian is 

Ho = -E E EffiWC^ + H - c -) ^EE n ia , 



i,j a,b—l (J 



i a—1 



(i) 



where c[ aa (ci acr ) is the creation (annihilation) opera- 
tor that creates (annihilates) an electron in ti g orbital 



a(= 1,2,3) with spin <r(=t, I) at site i, h ia = J2 a ^iao = 
E CT cJ oa .Ci acr , and H.c. means the Hermitian conjugate. 
Here, t 1 ^ (<j>) and fi denote the values of in-plane hopping 
integrals at an angle </> and the chemical potential, re- 
spectively. \i is determined so as to satisfy n c = A with 
n e being the total occupation number for the Ru ti g or- 
bitals. For simplicity, we neglect the effect of the spin- 
orbit interaction on the electronic structures; the validity 
is addressed in Sec. IV. 

The hopping integrals for Si - 2Ru04 with = 0° are 
determined^! so as to reproduce the FSs obtained in the 
de Haas- van Alphen effect^ the dispersions are given by 

en(fe, 0°) = — 1t\ cos k x — 2^2 n y — p, 
622 (fc, 0°) = —2^2 cos k x — 2t\ cos k y — fi, 
633(fc, 0°) = — 2t3(cos k x + cos k y ) — Ati cc 



£ab(fe,0°) 



%i cos ky — fi, (2) 

(3) 

k x + cos k y ) — At 4 cos k x cos k y — ji, 

(4) 

for a^i), (5) 



where tx/tz = 1.5, t 2 /t 3 = 0.2, and t 4 /t 3 = 0.4. In this 
tight-binding model, we neglect both the weak hybridiza- 
tion between the d xz and d yz orbitals and the difference 
of the CEF energy between the d xz / yz and d xy orbitals 
for simplicity^ In this work, we set t^ = 0.45 eV to make 
Wtot about 4 eV. Figures fJHl] show the band structure of 
the Ru tig orbitals, the FSs, and the DOS, respectively. 
In the next section, we describe the dependence of the 
dispersions on <fi for 0.5 < x < 1.5. 
The interacting Hamiltonian is 



Hmt = C/EE n iat n 



i a>b 



JhEE(2* 



ia ib 



) n ia n ib) 



(0) 



where U, U', and Jh are the intraorbital Coulomb 
interaction, the interorbital Coulomb interaction, and 
the Hund's rule coupling; ii a is defined as Sj a = 
(V 2 ) Ecr.o-' cloo-O'o-.o-'tW' with <r a>a , being the Pauli ma- 
trices. In the interacting Hamiltonian, we have neglected 
both the pair hopping J' and the transverse components 
of Jh for simplicity; the roles of these terms are discussed 
in Sec. IV. We thus consider the following Hamiltonian 
in the absence of the rotation of RuOg octahedra: 



int • 



(7) 



A. Effects of the Ca substitution 
on the electronic structures in 0.5 < x < 2 

As described in Sec. I, the Ca substitution induces 
the rotation of RuOe octahedra along the c axis in 
0.5 < x < 1.5. In order to take account of effects of 
the Ca substitution on the electronic structure, we as- 
sume that the main effects of the Ca substitution are the 
changes of the dp hybridizations due to this rotation. As 
shown in Fig. Q] (b), this rotation changes both the Ru- 
Ru lengths and the overlap integrals between the Ru Ad 
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FIG. 5: (Color online) Schematic picture of the in-plane dpir 
hybridization of the Ru d xz orbital to the O p z orbital in the 
presence of the rotation of RuOe octahedra. x and y (x' and 
y') are the coordinates in non-rotated (rotated) frame. The 
difference of a color in the d xz orbital represents that of the 
sign of the wave function for the d xz orbital. 
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FIG. 6: FSs for the effective model of x = 0.5. The dashed 
line represents the folded Brillouin zone in the presence of the 
rotation of RuOe octahedra. 



and O 2p orbitals keeping the Ru-0 bond lengths con- 
stant. The former leads to negligible effects on the elec- 
tronic structures since the angle dependence of the Ru- 
Ru lengths is cos 0, and the latter results in the change 
of the dp hybridizations linear in 0. 

We first explain the method to construct the tight- 
binding model by taking into account the effects of the 
rotation of RuOg octahedra on the electronic structures. 
As shown in Fig. [5J we can represent the wave functions 
for the Ru Ad and O 2p orbitals in the non-rotated frame 
(x, y, z) by those in the rotated frame (a/, y' , z'): 

-01 (r) = -0^ z (r) 

= ^d xU ,(r) cos <f>- ip dy , z ,{r)sm(j), (8) 

^a(r) = il>d VM ( r ) 

= tpd x , z , (r) sin + V> dy , 2 , (r) cos 0, (9) 
i>a(r) = ?pd xv (r) 

= ^d x , v , (r) cos 20 + tpd a , 2 _ yl2 (r) sin 20, (10) 
Mr) = ^d x2 _ y2 (r) 

= ~ i>d x , „, (r) sin 20 + if> d (r) cos 20, (11) 
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FIG. 7: (Color online) DOS for the effective model of x = 0.5. 
The dashed black line represents the chemical potential. 
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FIG. 8: (Color online) DOS for the special model. The dashed 
black line represents the chemical potential. 



and 



( r ) = ipp a , ( r ) cos - ip Py , (r) sin 0, (12) 

*Ppy ( r ) = M> Sm ^ + ^Py> ( r ) C0S <t>> ( 13 ) 

iM*-) = iM»0- ( 14 ) 



The dp hybridizations in the rotated frame are equal 
to those at = 0° since the Ru-0 bond length is un- 
changed^ we have Vd , , p , (0) = V\ Pz (0°) and so on. The 
finite dp hybridizations at = 0° are Vi Pz (0°), V 3p (0°), 
and V ipx (0°) along the x direction and V 2pz (0°), V 3p l (0°), 
and V\ Py (0° ) along the y direction, respectively. There- 
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fore, the dp hybridizations at are given by 

^W<« = ^(O°)cos0, (15) 
V 2p M = F lpj8 (O o )sin0, (16) 
Va P * (0) = ^ 3Ph (0°) sin cos 20 + V 4px (0°) cos sin 20, 

(17) 

^3 Ph (0) = U3pjO o )cos0cos20 + U 4 p x (O o )sm0sin20, 

(18) 

in the x direction, and 

^(0) =-V^(O°)sin0, (19) 
Vap.fa) = ^(O o )cos0, (20) 
V3 Px (0) = V£ Pa> (O°)cos0cos20 — V4 Pv (0°) sin sin 20, 

(21) 

V 3pv (0) = V 3px (0°) sin cos 20 + F 4pB (0°) cos sin 20, 

(22) 

in the y direction, respectively. Using the second- 
order perturbation processes with respect to these dp 
hybridizations, we obtain the following nearest-neighbor 
hopping integrals among the Ru t 2g orbitals: 



t^i B;x (0) = h cos 2 0, 



' n 

4^(0) 



i 



-t\ sin 20, 



t 



AB;x 
21 



ist>) 



AB;x 
12 



WO, 



t^' x ((j>) = t 2 — t\ sin 2 0, 
tf^' x ((j>) = t 3 cos 3 20 — tg cos 20 sin 2 20 
+ 2^6 cos 20 sin 2 20, 

in the x direction, and 

tf 1 B;y (0) = t 2 -Min 2 0, 



ti2 ;V (<t>) = - T^i sin 20, 



1 



AB-y, 



(0) = ti COS 0, 

t^' v (4>) = t 3 cos 3 20 — ^5 cos 20 sin 2 20 



2^6 cos 20 sin 20, 



(23) 

(24) 

(25) 
(26) 

(27) 

(28) 

(29) 

(30) 
(31) 

(32) 



in the y direction, respectively. Here, the superscripts 
A and B denote the sublattices, and t\, t 2 , t 3 , and 
are defined in Eqs. ©-flU- In deriving Eqs. (|23|) - (|32| . 
we have assumed that t\ in Eqs. ([2]) and ([3]) arises from 
the second-order perturbation process with respect to the 
dp hybridization of the Ru d xz / yz orbital to the O 2p z 
orbital: t\ is given by 



ti = 



^(0°) 2 



^(0°) 2 



E^O*) - E Pz (0°) E 2 {Q°)-E p ^°Y 



(33) 



where E a ((f)) and £'b(0) are the CEF energies at for 
the Ru t 2g orbital a and the O 2p orbital b, respectively. 



For the case with = 0°, we have V 3px (0°) = V 3py (0°) 
and E Px (0°) = E Py (Q°). In the presence of the rotation 
of RuC>6 octahedra, E Px (0) and E Py (0) can be different 
in principle; however, we have assumed E Pa . (0) = E Py (0) 
in deriving Eqs. (f23" |) - (j32"|) for simplicity. Similarly, t 3 in 
Eq. (0} is given by 



o *2 



^3p H (0°) 



v 3p Ao°) 



o\ 2 



E 3 (0°) - E Py (0°) E 3 {0°) -E Px (0°)' 



(34) 



In contrast, t 2 in Eqs. ([2]) and ([3]) arises from the direct 
hopping process between the d xz / yz orbitals since there 
are no hybridizations with the O 2p orbitals in between; 
we have assumed that the rotation of RuC>6 octahedra 
does not affect t 2 in Eqs. (J26J) and (|28j). In Eqs. flST} and 
(|32[) , is and te are the rotation-induced hopping integrals 
which arc defined as 



v iP Ao°) 



o >2 



^4 Pa (0°) 



o >2 



E 3 (0°)-E P M 
V 3Pv (0 o )V iPx (0°) 



E 3 (0°)-E Py &)' 

v 3p Ao°)v iPv (o°) 



E 3 {0°)-E P M E 3 (0°) - £ Ph (0) 



(35) 
(36) 



respectively. The former arises from the hybridizations 
of the d x '2. y >2 orbital in Eq. (jlpp with the p x i and p y i 
orbitals in Eqs. (fT2|) and (|T3|) . and the latter arises from 
the hybridizations of the d x 'y' and orbitals in Eq. 

(JTOJ) with the p y i and p x > orbitals in Eqs. (|T2|) and ([13")) . 

Using these hopping integrals Eqs. (j23|) — (|32 ]l . we ob- 
tain the noninteracting Hamiltonian in the presence of 
the rotation of RuC>6 octahedra as follows: 



H 



k 



6=1 s,s'=A,B 



E^'( fc > 



Wksa^ks'ba, (37) 



where 
e^(fc,, 

4 A (*v 
4 s (fc,< 



c ab 

-12 



(fc,0 



e^(fc,, 

4B 



AS 
£ 33 



AA 



BA 



(fe,, 



(fe,, 



(fe,, 
(fe,, 



4A 



(fe,, 



2/ 



2t 



AB ;x i 
12 



-12 



(fe. 



(38) 
(39) 

(40) 

((f)) cos k x — 2tf^' v ((f)) cos k y , (41) 

(42) 



4£4 cos k x cos k y — 
- 2tf^' x ((f>) cos k x 



AB;y 
11 



(0) cos k y , 



2t^ 2 B;x ((j)) cos k x - 2t^' ,y {4>) cos k y , 



AB;y, 



2t^' x (<fi) cos k x 



2t 33 i ' S '(0) COS fej,, 



AB 



(fc,0) =0 

(fc,-0). 



otherwise, 



(43) 

(44) 
(45) 
(46) 



The prime in the summation with respect to fc repre- 
sents the restriction within the folded Brillouin zone in 
the presence of the rotation of RuOg octahedra. For sim- 
plicity, we have neglected the effect of the rotation of 
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RuOg octahcdra on the next-nearest-neighbor hopping 
integrals for the Ru Ad orbitals. 

The rotation of RuOg octahcdra also induces the hy- 
bridization of the d xy orbital to the d x 2_ y 2 orbital al- 
though it is not included in Eqs. p3| - ([35]) . Since the 
energy level for the d x i_ y i orbital is higher than that for 
the d xy orbital, this hybridization leads to a downward 
shift of the d xy orbital. This effect is approximately taken 
into account as the difference of the CEF energy, A t2s , 



between the d xz / yz and d xy orbitals: 



HCEF = 



(47) 



In this work, we use At 2g as a parameter. Combining 
this term and the noninteracting and interacting Hamil- 
tonians, wc consider the following Hamiltonian in the 
presence of the rotation of RuOg octahcdra: 



int. 



(48) 



In order to discuss the electronic state for x = 0.5, we 
set (j) = 15° and (t 5 ,t 6 ,A t2g ) = (0,0,0.42) (eV) so as to re- 
produce the FSs obtained in the ARPES measurement^ 
Figures [6] and [7] represent the FSs and the DOS, respec- 
tively. In Fig. [6j the absence of the FS around k = (ir, 0) 
corresponds to the change of the FS for the d xy orbital 
from an electron pocket to a hole pocket. There are three 
main changes of the electronic structures due to the rota- 
tion of RuOg octahcdra: the reduction of the bandwidth 
for the d xy orbital, the increase of the pDOS for the d xy 
orbital near the Fermi level, and the change of the FS for 
the d xy orbital from an electron pocket to a hole pocket. 
The bandwidth and pDOS for the d xz / yz orbital are lit- 
tle affected. These results are qualitatively consistent 
with those obtained in the density-functional calculation 
within the LDA— or local-spin-density approximation.— 
We thus think that the present model is reasonable for 
the study of the electronic states in C^-aSr^RuC^ al- 
though the experimental value of <f> is about 12.8° at 
x = 0.54 

In addition, we study a special model in which the vHs 
for the d xy orbital is located on the Fermi level in order 
to clarify the role of the vHs in enhancing the effective 
mass. We use the dispersions in Eqs. (|38 ]) ~ (|46| setting 
4> = 15° and (t 5 ,t 6 ,At 2g ) = (0,0,0.3) (eV); the difference 
between this case and the case of x = 0.5 is that in the 
values of A t2 . Figure [8] shows the DOS in this case. 



B. Gutzwiller approximation 

In this section, we explain the GA for a PM state con- 
sisting of the degenerate d xz and d yz orbitals and the 
d xy orbital. Let us define the Gutzwiller-type variational 
wave function as 



(49) 



where l^o) is the ground state for Hq, Pj-i are the projec- 
tion operators, and gi are the corresponding variational 
parameters for the Ith configuration. There are 4 3 = 64 
configurations for each site j since we have three orbitals 
and two kinds of spin degrees of freedom. The variational 
energy is given by 



E gs = 



(tt G |ff|tt G ) 

{g7i (*g|*G> 



(50) 



where min{ 9( j denotes the optimization with respect to 
{gi}, and H is the total Hamiltonian given by Eq. (|48p. 

To calculate the expectation values, we use the 
GA22r— in which the dependence of the Slater deter- 
minants on the configurations is neglected. For example, 
the square of the Slater determinant, [det [7i-|-({r^. 1T })] 2 , 
for the electrons in a = 1 with a = j" can be approximated 
as 

[det C/ lt ({rf-})] 2 G 4 (n? t )*t(i - n?^"^ (51) 

where {»*^ 1T } is a set of sites occupied by these electrons, 
L is the number of sites, Ni-f is the number of these elec- 
trons, and is defined as an expectation value without 
the projection operators: 

»?t = T E<*o|%i t |*o) = 7 EflWo- (52) 

Using the GA, we can rewrite the denominator in Eq. 
([5D]1 as 

<*g|* g > = e e [n. 9 r i({ ^ }) ] 
x ft n[dctt/ Q ,(KV})] 2 



^ A e E(n^ 2r! ) 

{iWfr,} i 



(II iv 



■P(L;{N ar7 }), 

(53) 



with 



P(L;{N aa }) = l J(n° a<T ) N ^(l-n° aa ) 



\L — N a 



(54) 



Here, E{./v } denotes the summation with respect to 
all the possible values N aa (a = 1,2,3, a = ^,1), Ti is 
the number of sites with the Ith configuration, and the 
summation E{T(} ^ s taken over all the possible values Ti 
(I = — 63) under the constraints 



j 1=0 



63 

E r ; = l, 

1=0 

63 

E E r ; n /acr = J2 N °-<7, 

1=0(1=1,2 a=l,2 
63 

E Fini 3cr = N 3a , 

1=0 



(55) 
(56) 
(57) 
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where ni aa is the number of electrons in orbital a with 
spin a for the Zth configuration; as we shall show, n\ aiJ 
plays important roles in calculating the expectation val- 
ues within the GA. 

Since the occupation number for each orbital is same in 
A and B sublattices, we do not need to take care of the 
sub-lattices in the GA. In other words, the formulation 
for a PM state consisting of the t 2g orbitals without the 
sublatticc indices is sufficient to calculate the expectation 
values for a PM state consisting of the t 2g orbitals with 
the sublattice indices. 

(iV^ + N 2a ) and A^3 CT are conserved quantum numbers 
although Niu and N 2a are not conserved quantum num- 
bers due to the finite hopping integral between the d xz 
and d yz orbitals. In general, the values of (Ni a + N 2<y ) 
and N 3a may change in the presence of interactions for a 
system with orbital degrees of freedom. In the following, 
however, we assume that these values are equal to those 
in the absence of the interactions; e.g., for the effective 
model of x = 2, {N Xa + N 2lJ )/L = (n? + n9.) = 2.66, 
N 3lJ /L = n§ = 1.33. In Sec. IV, we address the effect of 
the variation of the occupation numbers due to electron 
correlation. 

In the thermodynamic limit, we can approximate the 
summation with respect to Ti in Eqs. (fS"3"f to their largest 
terms: 

L\ 



and 



GA 



gi*g; 



(n 



.21', 



9i 



aim 

i 



P(L;{N aa }), (58) 



where {T;} is the set of {r;} that gives the largest term. 
Ti is given by 



Ti = Lgfcx.p(l + A + EE Ko-niaa), 



(59) 



where (Ao, {A acr }) are the Lagrange multipliers^ deter- 
mined so as to satisfy the constraints Eqs. (f55|) -([57 |) for 
{f For simplicity, we assume that the following rela- 
tions hold in the thermodynamic limit: 



(60) 
(61) 



N at _ 




L 


L 


N la 


N 2a 


L 


L 



Correspondingly, the Lagrange multipliers satisfy 
A a t = A a j,, 

A2cr = Al CT . 

Similarly, the expectation values of ni aa -niba 
ct atT Cjba can be rewritten within the GA as 

l*G> 

^» ^2 ^ 9u. (Pi;h fiia^n ibcr > Pj-ijo (j] 9^' 
U {r;> i 



(62) 
(63) 

and 



<nr;- 



j±P(L-l;{N' aa }), 



(64) 



GA — " 



x (n3ro^rffT^-2;ra), 



(65) 



with the constraints 



63 



E r{ = l - 1, 

1=0 

63 

E E r> a(7 = E Ka 

Z=0a=l,2 a=l,2 



= E (Naa ~ n hai j), 
a=l,2 



63 



E I> 3 „ = 
1=0 



N 3 a - n h3o 



and 



63 

E rf = l - 2, 

1=0 



(66) 



(67) 



(68) 



(69) 



63 



E E lYrw 

1=0 a=l,2 



63 

E ^>l3o 
1=0 



: y n" 

a<7 

a=l,2 

: E (Nacr ~ n h aa ~ n^aa) 
a=l,2 

: E (Naa ~ ni'. aa - n V , aa ), 
a=l,2 3 



n'L 



N 3a - ni.3a - ni j3a 

N 3a - 1lE3a - n>l'3a- 



(70) 



(71) 



Here, the summation with respect to li in Eq. (|64[) is 
taken for the configurations of site i in which there are at 
least an electron in orbital a with spin a and an electron 
in orbital b with spin a' . In Eq. (|65j) . li and lj are the 
configurations of sites i and j after the hopping process 



i aa Cjba in which there are no electrons in orbital b with 



spin a at site j and there is at least an electron in orbital 
a with spin a at site i; l\ and are similar configurations 
of sites i and j before the hopping process. The prime 
and double prime in the summations of Eqs. (|64[) and 
(p5|) represent these restrictions. 

In the thermodynamic limit, we approximate the sum- 
mation with respect to FJ in Eq. (|64p to their largest 
terms for each configuration of li and the summation 
with respect to T'/ in Eq. (|65jl to their largest terms 
for the configurations of li, lj, l' it and lj. Thus, these 
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expectation values become 

(*g| l*G> 

^ 9^ (Pi-h fliaaflibrr' A;i< )o (T\_ 9, 
h I 

„ (i-l)' 



GA 



2i : 



01 f 



-^P(L-l;{iV^}), 



(72) 



and 



(*g|cLcj6<t|*g) 



GA 



XI 91^1,91'^ (hhPjlhclaa^baPi-l'.Pjd^O 



(Tin 

i 

where T[ and r" are given by 



(73) 



and 



(*G|cLc j&CT |^ G ) 



GA 



(*g|*g) 

1 ^-V' » 4. » - 

^ 2^ aughgi'.gi'jiPi-MPj^Ciaa^Pi^Pj^io 



v- 



li,U,E,l'. 



Tl9\ 



^(Tim P(i _ 2;WJ) 



air;'!) p{L-{N aa }) 
i 



(85) 



where AT J and Arj' are defined as Arj = T' l -T l and 
AT;' = T'{ - respectively. From Eqs. (^HEZI) for 
{I>}, Eqs. (|66])-([68l) for {rj}, and Eqs. dp-dTT]) for 
{r"}. we have the following constraints for Arj and Ar": 



f i = Lgf exp(l + K,nia.) , (74) 



a — 1 a 
3 



F" = Lgf exp(l + Aq + E E Ka n ia 

a—1 <7 



(75) 



respectively. The Lagrange multipliers, (Aq, {A^}) and 
(Aq, {A"^}), are determined so as to satisfy the con- 
straints Eqs. (jMl-dSll) for {fj} and Eqs. (JeHl) dZU for 
{r"}, respectively. For simplicity, we also assume that 
the following relations hold in the thermodynamic limit: 



and 



L 
L 



N" 

"at 



N'L 



N' 

IT' 

L ' 



n'L 



L L 

Correspondingly, the Lagrange multipliers satisfy 



A„* — A„ 



A^x — A 1(T , 



(76) 
(77) 

(78) 
(79) 

(80) 
(81) 



and 



A 2<T 



(82) 

C (83) 
Using Eqs. (|58|) . (f72|) . and (|73)l . we obtain the normal- 
ized expectation values of fiiacrfiiba' and c\ aa Cjba within 
the GA in the thermodynamic limit: 

(* ' G\hia<jhib<ji\^ g) 



GA 1 
-> - 

L 



(*g|*g) 

gh( P i\hf l ia<jnibcr'Pi:l i )o 

V 7 



2Ar 



P(L-l;{j\£J) 



(nrj!) P(i;{iV aCT }) 



(84) 



63 

63 

E E At> a(7 

;=0 a=l,2 
63 

E Af \n l3a 

1=0 



E n haa, 
a=l,2 



and 



63 

E Arf 

63 

E E Af;'n ;aCT 

i=0 a=l,2 



63 

E Af >/n l3o 

1=0 



= -2, 



(86) 
(87) 



E ( n haa + ni. a<J ) 
a=l,2 

E ( n l'.aa + ni'.ac), 
a=l,2 1 3 

("/iScr + ni j3a ) 

(ni'. 3 a + ni'. 3a ). 

P(L-l;{N' aa }) ■ 



(89) 



(90) 



(91) 



Calculating (Pi-uhiaahibe' Pi-i^o and — ^ 



{iv OJ }) inE q- 

P(L-2;«;}) 



and l):iJ):i J ' : l, T r ii ,,J' l: i : l'j:r i o and P(L . {jW) 
in Eq. (|85[) explicitly for each possible configurations, we 
find 



(^Glniacrfliba'l^G) GA 1 



(*g|*g) 



2Ar;^ ; 
3z 



OK, 



(nr;o' 



(92) 



and 



9 



(*g|cT oct c j - 6ct |* g ) ga 1 



(*g|*c 



2Ar"\ / 



(nr, 



li,U,l'.,l'. ' 

In addition, we have the following relation in the thermodynamic limit: 

(fi)l 



(r ; + Ar,)! 

Using this relation, we can rewrite Eqs. (|92p and as 



(*G|«<acr«iI><r'|*G) GA 1 



L 



E'^n 



(*G|cLcjto|*G) GA 1 



(*g|*c 



L 2 



E" 



(Tj) 



-Ar! 



-Ar, 



1 



iaa 



n 



(i-»2K Wtf 



-Ar' 



(93) 

(94) 

(95) 
(96) 



respectively. Substituting Eq. (|59p into these equations and using relations L Ej Ar < = L and L Ar ' = L 2 from 
Eqs. (J5BJ) and ASS), we find 



|*g) ga 



^'.g ; 2 i exp[-^(l + A )Af;- £ WtwAT; 



E -^CXP (l + ^0)+ E 



E 



(97) 



(*Gl4 g g jtol*G) 

(*g|*g) 



GA 



E 91^91'^ {ct a *Cjba)o . _ Kr-p exp[- ^(1 + A )Af;'] cxp[- ^ A ' ff »Tii v»Arf 

h,U.l'..V. \ a) b i I. a', a" 



- £ 

U.l^l'..V. 



9h9l i 9V i 9V j (CiavCjbalo-r "oTTcr ex P 



2(1 + A ) 

eXP [ 2 E ^o,'a'{lM ia 'a' + n lja , al + n Via , a , + ni'.a'a') 



l ' 3' j 

In deriving Eq. (|98|) . we have also used the following relation due to Eq. (|63)) : 

E A aCT (n ; / 

acr T TT'l'.aa) ~ E A aff (n z . aCT + n iaCT ). 
a=l,2 3 o=l,2 

Finally, the variational energy becomes 

2 



„ GA . 

h ss — > mm 



-E E E CbW (1 _ n0)n o (4 g g j^)o E 

i,j o,6=l<r=t4 V 1/1 U,U,l'.,l'. 

J t J 



V- 



E E E ffi W (i_ n o )n o (4^)o E 



„ wr/.r z T,/.r ; 



3 ' i' 3 



L 2 



(98) 



(99) 



^EEEE £'% + ^'EEEE'%- ^EEEE'I 



i a— b a="X cr'— J, l± 



i a>b a— a' l± 



(100) 
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TABLE I: The configurations, the wave functions, the values of Hint, and the variational parameters gx k for a PM state 
consisting of the degenerate d xz and d yz orbitals and the d xy orbital. and t (4.) in the wave functions mean there are no 
electrons in orbital a with spin a and there is a spin-up (spin-down) electron in orbital a, respectively. 



Wave functions |1 t, 1 4-! 2 t> 2 4-! 3 t> 3 4-) value of H int gx k 



X 








|0, 0; 0, 


0;0,0) 











9x 


Xi 


It.c 


;0,0 


0,0), 


|0, 4-; 0,0; 0,0), 


|0,0;t, 0; 0,0) 


10, c 


;0,4 


;0,0) 





9x 1 


x 2 








|0,0;0,0;t,0), 


|0,0;0,0;0,t) 











9x 2 


x 3 








|t, I; 0,0; 0,0), 


|0,0;t,l; 0,0) 








u 


9x 3 


x 4 








|0,0;0, 


0;t,t) 








I! 


9x 4 


x B 








|t, 0; 0,|; 0,0), 


|0,|;t,0;0,0), 








V 


9x 5 


X 6 


IT, C 


;0,0 


0,1), 


|0,4-;0,0;t,0), 


|0,0;t,0;0,t) 


10, c 


;0,4 


;t,o) 


U' 


gx 6 


X- 








|t,0;t,0;0,0), 


|0,4-; 0,4-; 0,0) 








U' - J H 


9x 7 


X s 


It.c 


;0,0 


t,o), 


|0, 4-; 0,0; 0,4-), 


|0,0;t, 0;t,0) 


10, c 


;0,4 


;0,4-) 


U' - J H 


9x 8 


Xg 


It, 4 


-;t,o 


0,0), 


|t, 4-; 0,|; 0,0), 


lt,0;t,l;0,0) 


lo, 4 


-1 T? ■> 


;0,o) 


U + 2(7' - Jh 


gx 9 


X10 


It, 4 


-;0,o 


t,o), 


It, 4-; 0,0; 0,1), 


|0,0;t,l;t,0) 


10, c 


;t,4 


;0,|) 


U + 2U' - J H 


gx 10 


Xn 


It.c 


;0,0 


t,4->, 


|0,|;0,0;t,4-), 


|0,0;t, 0;t,t) 


10, c 


;0,4 


;t,4-) 


U + 2(7' - J H 


to. 


X12 


|t, 0; 04; 0,4.), |0,4 


-;t,o 


t,o), 


lt,0;0,|;t,0), 


|0,t;t,0;0,t) 


It.c 


;t,C 


;0,|), 


|0, 1; 0, |; t, 0) 317' -Jh 


3^12 


X13 








lt,0;t,0;t,0), 


|0,4-; 0,4.; 0,4,) 








3(7' - 3J H 


gx 13 


X\4. 


It, 4 


;t, 


0,1), 


lt,t;0,i;t, 0), 


It, 0: t, 4.; 0, 1) 


10, J 


;t,4 


;t, 0) 


U + 5(7' - 2 J H 


14 


X15 








lt,0;0,|;t,l), 


|0,t;t,0;t,t) 








U + 5(7' - 2 J H 


gx ls 


XlQ 


It, 4 


-;t,o 


t,o>, 


It, t; 0,|; 0,|), 


|t,0;t.4-;t,o) 


10, 4 




;0,|) 


U + 5(7' - 3 J H 


gx 16 


Xl7 








lt,0;t,0;t,l), 


|0,t;0,|;t,t) 








U + 5(7' - 3 J H 


gx 17 


Xis 








|t,t;t, 


4-; 0,0) 








2(7 + 4(7' - 2J H 


gx ls 


Xig 








lt,t;0,0;t,l), 


|0,0;t,l;t,t) 








2(7 + 4(7' - 2J H 


gx 19 


X20 


It.c 


;t,4- 


t,4->, 


|0,t;t,l;t,l), 


lt,t;t, 0;t,t) 


It, 4 


-;0,4 


;t,4) 


2(7 + 8(7' - 4J H 


gx 20 


X21 








lt,t;t,l;t,o), 


It, 4-; t. 4-; 0,4-) 








2(7 + 8(7' - 4J H 


gx 2 i 


X22 








|t,4-;t, 


4-;t,t) 








3(7 + 12(7' - 6 J H 


gx 22 



Here, the restrictions about the summations with respect to Zj, lj, and in the first and second terms are those 

for the hopping process c\ aa Cjba (a, 6 = 1,2) and hopping process ct 3o .Cj3 CT , and the restrictions about the summations 
with respect to Zj in the third, fourth, and final terms are those for the interaction terms U, U', and Jh- In deriving 
Eq. (|100j) , we have neglected the term of Hqef in the variational energy since this term only gives a constant energy 
shift which is independent on the strength of the onsite interactions within the present treatment. Note that E gs 
can be written as the products of the component for orbital a with spin a since we have neglected both J' and the 
transverse components of Jh in the interacting Hamiltonian as denoted in Sec. II. 

For simplicity of the numerical calculation, we classify the possible 64 configurations into 23 groups in which the 
energies are same. These are shown in Table HI where Xk denotes the optimized number of sites with the configuration 
in the k group; e.g., X\ = Y\ for I = 1 — 4, and so on. Denoting Xk = X^/L, we obtain 



GA 



mm 

X 3 ,X4,--- ,X 22 



~E E E a a&(x 3 ,X 4 ,--- ,X22)t l a Z((t>)(c\a*Cjb*)o 
i.j a, 6—1 fx 

+ (7(2x3 + x 4 + 4x9 + 4xio + 4.Tn + 4£i4 + 2xi 5 + ix 16 + 2xyj + 2x 18 + Ax 19 + 8x 2 o 
+ 4.T21 + 3x22) 

+ U' (2x 5 + 4x 6 + 4x 9 + 4xio + 4xn + 12x i2 + 12x i4 + 6x i5 + 8x i6 + 4xi 7 + 2xi 8 + 4xi 9 

+ 16x 20 + 8x21 + 6x 22 ) 
+ (V - J H )(2x 7 + 4x 8 + 4x 9 + 4xio + 4xn + 6xi 2 + 6x13 + 8x14 + 4xi 5 + 12xi 6 + 6xi 7 

+ 2xi8 + 4xi9 + 16x20 + 8x21 + 6x22) • (101) 



11 



Here, q a b(x3, %4, • • • , £22) is the RF of the kinetic energy for the Ru t 2g orbital, which satisfies 

qu(x3,XA, ■■■ ,£22) = 912(^3, xa, •• ■ ,£22) = 921 (^3, £4, • • • ,£22) = 922(^3, x±, ■ ■ ■ ,£22)- (102) 
The RFs for the d xz / yz and d xy orbitals are given by 



'In 



i(ar 3 ,a;4, • • • ,x 22 ) = , , [v^i(v^o + v^3 + ^+V^7 ) + V^eiV^ + a/^io + 2^/x^ 
n i(t n 1 ) 



and 



X8(y/X2 + y / Xw + \fxyi + Vxi3 ) + \fx$(y/xi + \fx~l + y/x~7 + y/xi8 ) 
+ \fi^(y/X4, + V^15 + V^VT + \/^19 ) + V^uiV^O + 2 \/^12 + y/xzi ) 

+ v/^ie(\Aio + \Ai2 + V^ia + \/^2i ) + v / ^(v / ^i5 + >Air + V%i9 + \/^22 )] 2 , 

(103) 



933(^3, X 4 , • • • ,X 22 ) 



»§(l-n°) 



qT [v^Cv^O + ) + 2y/x~[(y/x^ + ) + 2 y /x~^{ y /x 3 ~ + V^19 ) 



+ 2y/xTl{\/x^ + ,/x^ ) + 2y/x^(y/x^ + y/xTE ) + \[x^{ypX^2 + V^13 ) 

+ v / ^7(V^i2 + ) + 2^/a^( % /iH + v^ie ) + 2^fx^{^xTl + ^/x^ ) 

)] , (104) 

respectively. Note that the optimization with respect to is equivalent that with respect to {xk} due to the 

relation Xk = g\ k exp(l + Ao + ^2 a a \ aa ni aa ). In the numerical calculations, we use the following constraints instead 
of determining the Lagrange multipliers Ao and A acr : 

1 = x + 4ei + 2x 2 + 2x 3 + X4 + 2x 5 + 4x 6 + 2x 7 + 4x$ + 4x 9 + 4x 10 + ixn + 6x 12 + 2x 13 + 4xi4 + 2x 15 

+ 4.ti 6 + 2xn + xis + 2x w + 4x 2 o + 2x 2 i + £22, (105) 
= x\ + xz + x 5 + x 6 + x 7 + xs + 3xg + 2x w + xn + 3xi 2 + £13 + 3xi 4 + x 15 + 3x ie + x 17 + xi$ + xig 

+ 3x20 + 2x 2 i + X22, (106) 
n 3 = x 2 + X4 + 2x e + 2x s + 2x w + Axu + 3x 12 + X13 + 2xu + 2x 15 + 2xi 6 + 2x 17 + 2x lg + 4x 2 o + £21 + ^22- (107) 



III. RESULTS 

In this section, we show the numerical results of the GA 
for three cases with the effective models of x — 2 and 0.5 
and the special model. The variational energy Eq. (jlOip 
is numerically minimized with respect to {xk} under the 
constraints (|105[) - (|107l) by Powell's method^ which is 
one of the numerical methods to minimize a function with 
more than one variable. In this work, we use the value 
of U as the parameter, and set U' = U — 2 Jh, Jh = U/A, 
and Wtot =4.1 eV, which is obtained in the effective 
model of x = 2. 



A. Mass enhancement for the effective models of 
x = 2 and 0.5 

We first show the results for the effective models of 
x = 2 and 0.5. Figures l9l and flOl represent the RFs of 
the kinetic energy for the Ru i 2 g orbitals as a function 
of U /Wtot for these models. We see from Fig. [H] that the 
RFs for the d xz / yz and d xy orbitals are nearly the same 



for the case of x = 2. This is probably originated from the 
fact that these orbitals have nearly the same occupation 
numbers (i.e., n\ = nSj = n§ = 4/3) in the absence of 
the interactions. On the other hands, for the case of 
x = 0.5, we find from Fig. [TDIthat the RFs for the d xz / yz 
and d xy orbitals are different, and that the increase of 
U /Wtot leads to the large difference between these RFs. 
The occupation number for each Ru £ 2£ , orbital becomes 
(n?,ng,rag) = (1.17,1.17,1.66) for the case of x = 0.5. 
This change of the occupation numbers results mainly 
from the downward shift of the d xy orbital since we have 
{n\,n%nl) = (1.32, 1.32, 1.35), which arc little different 
from those for the case of x = 2, for the model setting </> = 
15° and A t = cV in Eq. ([S5]). By using the analogy 
with the result for the single-orbital Hubbard model, the 
difference of the RFs between the cases of x = 2 and 
x = 0.5 will be due to this change of the occupation 
numbers approaching the integer values towards x = 0.5, 
which is expected in the usual Mott transition; in the 
present case, the occupation numbers expected in the 
Mott insulator are 1 for the d xz / yz orbital and 2 for the 
d xy orbital, respectively. 

In order to compare our results with the cxperimcn- 
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0.2 0.4 0.6 0.8 1 1.2 1.4 

u/w tot 

FIG. 9: (Color online) RFs for the d xz i yz and d xy orbitals as 
a function of U/Wtot for the effective model of x = 2. We set 
Jh = (7/4. 
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FIG. 11: (Color online) Total RFs as a function of U/W tB 
for the effective models of x = 2 and 0.5. We set Jh = (7/4. 
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FIG. 10: (Color online) RFs for the d xz / yz and d xy orbitals as 
a function of U/Wtot for the effective model of x = 0.5. We 
set J H = (7/4. 
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FIG. 12: (Color online) RFs for the d xz / yz and d xy orbitals as 
a function of U/Wtot for the special model. We set Jh = (7/4. 



tally observed mass enhancement, we define a total RF, 



<ztot = y 3 E^L. (io8) 

which estimates the inverse of the mass enhancement. 
Figure [Til shows the total RFs as a function of U/Wtot 
for the effective models of x — 2 and 0.5. We see that g tot 
for the case of a; = 0.5 becomes much smaller than that 
for the case of x — 2 as U /Wtot increases. Therefore, our 
results suggest that moderately strong Coulomb interac- 
tion and the modifications of the electronic structures for 
the Ru t2 g orbitals due to the rotation of RuOg octahe- 
dra lead to mass enhancement for x = 0.5 than that for 
x = 2. 



B. Primary effect of the Ca substitution 



on mass enhancement around x = 0.5 

In order to clarify the role of the Ca substitution in 
forming the HF, we calculate the RF for the special 
model. Figure [T^] shows the RFs for each Ru t2 g orbital. 
We see that the RF for the d xz / yz orbital is slightly larger 
than that for the case of x = 0.5, while the RF for the d xy 
orbital is slightly smaller. The former results mainly from 
the decrease of the occupation number for the d xz / yz or- 
bital from 1.21 to 1.17, and the latter results mainly from 
the increase of the occupation number for the d xy orbital 
from 1.57 to 1.66; both changes of the occupation num- 
bers arise from the downward shift for the d xy orbital, 
which is induced by the rotation of Ru06 octahcdra. 

Figure [13] shows the total RFs for both the effective 
model of x = 0.5 and the special model. We find that 
the inverse of q tot for the effective model of x — 0.5 is 
larger than that for the special model, and that Coulomb 
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FIG. 13: (Color online) Total RFs as a function of U/Wtot 
for the effective model of x = 0.5 and the special model. We 
set J H = (7/4. 



interaction enhances this difference of the inverse of gtot ■ 
Therefore, our result indicates that the vHs for the d xy 
orbital plays a secondary role in enhancing the effective 
mass around x = 0.5; the primary one arises from the 
change of the occupation numbers approaching the in- 
teger values. Note that these occupation numbers are 
different from those expected in the OSMT but same to 
those expected in the usual Mott transition, where the 
occupation numbers for the d xz / yz and d xy orbitals are 1 
and 2, respectively. 



IV. DISCUSSION 

A. Comparison with previous theoretical studies 

We first remark on the physical meaning of the en- 
hancement of the effective mass obtained in the GA. For 
simplicity, we consider a PM state consisting of a single 
orbital. The following argument is based on the discus- 
sion in Ref. 33. One of the origins of the HF behavior is 
the criticality approaching the Mott transition. Namely, 
the formation of the HF results from separation of energy 
scale between spin and charge degrees of freedom due to 
electron correlation; on-site charge fluctuation suppresses 
at the higher scale corresponding to U, and the lower one, 
Tp, associates with the local spin fluctuation. Below Tp, 
a Fermi-liquid description applies; the mass enhancement 
is inversely proportional to a ratio of Tp to the noninter- 
acting Fermi temperature. Although the analysis using 
the GA restricts to the discussion about the properties 
for the ground states and can not address the possibil- 
ity of the above separation of energy scale, the GA can 
analyze the effective mass for systems with a different pa- 
rameter. A similar argument can apply to a system with 
orbital degrees of freedom. Therefore, we think that the 
analysis using the GA can capture the tendency of mass 
enhancement in Ca2_ a; Sr 2 ;Ru04 for 0.5 < x < 2. 



In the present calculation, we have not discussed Mott 
transition. From a theoretical point of view, the pre- 
vious DMFT calculation^ based on the models for 
Ca2- 2; Sr 2 ;Ru04 in 0.5 < x < 2 has shown that the criti- 
cal value of U for the Mott transition, U c , at x = 0.5 is 
larger than 6 eV (i.e., U/W to t > 1.4) at T = 0.02 eV. In 
principle, the Gutzwiller-type variational wave function 
gives the Mott transition as a function of U, when charge 
fluctuations, which play an important role in the vicinity 
of the Mott transition,— are included.— ~— However, the 
present study using the GA is inappropriate to address 
the possibility of the Mott transition. Thus, we have 
shown the results only for the case with U/Wtot < 1-4 
where the system will remain metallic. 

The present formalism of the GA neglects the varia- 
tion of the occupation number for each orbital due to 
electron correlation for simplicity. The previous DMFT 
calculation^ based on the models for Ca2- 2; Sr :E Ru04 
with 0.5 < x < 2 has observed a drastic variation of 
the occupation numbers (~ 20 — 30%) for values of U 
near U c , while the variation is about 10% (a few %) for 
U ~ 0.8t/ c (U ~ 0.7U C ). This indicates that the variation 
of the occupation numbers due to electron correlation is 
important near the Mott transition. In this paper, we 
have discussed the HF behavior only in the regions where 
the effects of the variation of the occupation numbers is 
small. 

Although it will be necessary to include the variation 
of the occupation number for each orbital for the quanti- 
tative argument, we think from the following arguments 
that the obtained tendency of the mass enhancement in 
0-5 < x < 2 docs not change qualitatively even if the 
variation of the occupation number for each orbital is 
included. According to the DMFT study*^ for the three- 
orbital Hubbard model, U c decreases when the occupa- 
tion number for each orbital in the absence of interactions 
approaches to an integer value. Simultaneously, the oc- 
cupation number for each orbital changes as a function of 
U/Wtot- Since the decrease of U c leads to the increase of 
the effective mass at a fixed value of U/Wtot and Jn/U, 
the mass enhancement will be largest at x = 0.5 where 
the occupation numbers in the absence of interactions 
are nearest to integer values in 0.5 < x < 2. Therefore, 
our results about the mass enhancement in 0.5 < x < 2 
will not change qualitatively even if the variation of the 
occupation number for each orbital is included. 

DMFT studies^r— proposed that the Hund's rule cou- 
pling plays an important role in stabilizing a metallic 
state with large effective mass. For example, a DMFT 
study^ for the degenerate three-orbital Hubbard model 
showed that the increase of Ja/U leads to a decrease of 
U c at i-filling and an increase of U c at i-filling. (Note 
that there are three electrons per a site at ^-filling and 
there is one electron per a site at ^-filling.) From this 
result, they claimed that the increase of Jn/U extends 
a region for the metallic state with large effective mass 
at the filling away from ^-filling; this metallic state is 
called Hund's metal. However, in Ca2_ a; Sr a ;Ru04, the 
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total filling number does not change in 0.5 < x < 2 (i.e., 
|-filling) and the value of Jh_/U does not change either. 
The main changes due to the Ca substitution are both 
the bandwidth for each Ru t2 g orbital and the occupa- 
tion number for each Ru t2 g orbital. Therefore, the effect 
of the Hund's rule coupling is not important for the dif- 
ference between x = 0.5 and 2. Instead, the criticality 
approaching the Mott transition plays a more important 
role in enhancing the effective mass in Ca2- £C Sr ;E Ru04 
towards x = 0.5. 

We now address the effect of the neglected terms in the 
interacting Hamiltonian (i.e., J' and the transverse com- 
ponents of Jh) on the electronic states. According to the 
previous study 2 ^ of the PM state for the two-orbital Hub- 
bard model in the GA, these terms slightly stabilize a PM 
metal against a nonmagnetic insulator. This result sug- 
gests that the effect of these neglected terms little affects 
the properties for the PM state studied in the present 
study within the GA. On the other hand, the transverse 
components of Jh will strongly affect the stability of the 
magnetically ordered states. Therefore, these neglected 
terms will not change the obtained tendency of the mass 
enhancement in PM states. 

According to the density-functional calculation for 
Ca 2 _ a; Sr a; Ru04 in 0.5 < x < 2 within the LDA^ the 
total bandwidth for the Ru t 2g orbitals becomes smaller 
for x = 0.5 than that for x = 2. Although this reduc- 
tion is partially included in our calculation through the 
4> dependence of the dispersions, it seems that these are 
other effects which reduce the total band width. Thus, 
if this effect is fully included, the difference between the 
values of g t ~J for x = 2 and 0.5 will become larger than 
for the present calculation. 

In our analysis, we have assumed that the roles of the 
O 2p orbitals are only to change the dp hybridizations. 
However, the previous theoretical work^ based on the 
dp model for Si~2Ru04 has proposed that Coulomb in- 
teraction for the O 2p orbitals plays an important role 
in stabilizing spin-triplet superconductivity. We expect 
that the Coulomb interaction for the O 2p orbitals leads 
to a larger mass enhancement than the present calcula- 
tion. The more detailed theoretical study about the role 
of the O 2p orbitals is deserved. 

There is a theoretical proposal that the vHs for the 
d xy orbital plays a primary role in forming HFs around 
x = 0.5 on the basis of the density- functional calcula- 
tion within the LDA4£ This work has proposed that the 
rotation-induced hybridization of the d xy orbital to the 
d x 2_ y 2 orbital causes the magnetic instability due to the 
nesting of the FS for the d xy orbital, and that the instabil- 
ity will lead to the mass enhancement around x = 0.5. In 
contrast, our results suggest that the primary role arises 
from the criticality approaching the usual Mott transi- 
tion resulting from the change of the occupation numbers 
for the Ru t 2g orbitals due to the downward shift of the 
d xy orbital, and not from the vHs. It will be necessary 
to study the effect of the vHs on the formation of HFs 
around x = 0.5 more systematically. 



B. Correspondence with experimental results 

We first discuss the roles of the Ru e g and O 2p orbitals 
in determining the electronic states for Ca2- a: Sr a; Ru04 
in 0.5 < x < 2. In this study, we have taken account of 
the effects of these orbitals as the changes of the dp hy- 
bridization and the CEF energy for the Ru t<i a orbitals; 
the latter is the downward shift of the d xy orbital due 
to the hybridization with the d x i_ y i orbital. There is 
no experimental evidence that the Ru e g orbitals play 
an important role in determining the electronic states 
except the possible change of the CEF energy through 
the hybridization of the Ru i 2g orbitals; thus, our treat- 
ment about the Ru e g orbitals will be valid. In con- 
trast, a polarized neutron diffraction measurement for 
Cai.5Sro.5Ru04 has observed a field-induced magnetic 
moment on the in-plane O ions, which is about 20% of 
that for Ru ions4£ This result indicates that not only the 
dp hybridizations, but also the Coulomb interaction for 
the O 2p orbitals will play non-negligible roles in deter- 
mining the electronic states. It is thus necessary to study 
the role of the O 2p orbitals systematically. 

Let us remark on the role of the spin-orbit interac- 
tion, which has been neglected in this work. The exper- 
imentally observed FS o 20 ' 32 for x = 2 and 0.5 are repro- 
ducible by the density-function calculations without the 
spin-orbit interaction i 27 ' 28 ' 33 Therefore, the spin-orbit in- 
teraction will be negligible in determining the electronic 
states. 

We have also neglected the effect of the disorder in- 
duced by the Ca substitution in this work. A measure- 
ment with a SQUID magnetometer for Cai.sSro.5Ru04 
has observed glassy behavior (i.e., the time-dependent 
magnetization)^ which is similar to the behavior ob- 
served in Cao.95Sro.o5Ru03.— This behavior is related 
to the disorder by the Ca substitution. It is thus nec- 
essary to include the effect of the disorder on the elec- 
tronic state in order to discuss the electronic states for 
Ca2- 2; Sr :Z ;Ru04. This remains as a future problem. 

Let us discuss the role of the vHs for the d xy or- 
bital. Experimentally, substitution of La 3+ for Sr 2+ in 
Sr 2 Ru04 leads to the downward shift of the vHs towards 
the Fermi level without any structural distortions^ In 
this case, the coefficient of the electronic specific heat for 
Sri. 8 La .2RuO4 reaches 1.3 times of that in Sr 2 Ru04. 5 ° 
This mass enhancement is mainly due to the increase of 
the DOS by the vHs for the d xy orbital; the density- 
functional calculation within the LDA and the ARPES 
measurement support this mechanismi 49 i 51 However, in 
the case of Ca 2 - a; Sr 2: Ru04 with x = 0.5, the ARPES 
measurement^ has shown that the vHs is located below 
the Fermi level. Therefore, there must be other factors 
other than the vHs in enhancing the effective mass in 
Ca2- 2; Sr 2 ;Ru04; our results indicate that one of the fac- 
tors is the criticality approaching the usual Mott transi- 
tion. 

We next remark on the roles of the spin fluctuations. 
Around x = 0.5, the resistivity behaves metallic (i.e., 
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dp/dT > 0), and the spin susceptibility shows Curie- 
Weiss behavior, as described in Sec. I. These experi- 
mental facts indicate that the mode-mode coupling for 
the spin fluctuations plays an important role in the elec- 
tronic states around x = 0.5 since the enhancement of the 
mode-mode coupling generally leads to the Curie- Weiss 
behavior^ The study taking account of the mode-mode 
coupling is a remaining future problem. 

There are several experimental results which indicate 
that ferromagnetic spin fluctuation also plays an impor- 
tant role in the electronic states around x = 0.5. One 
of the examples is the enhancement of the Wilson ra- 
tion towards x = 0.5Z This result indicates that the sys- 
tem is a nearly ferromagnetic metal near x = 0.5i 53 i 54 
In addition, the inelastic neutron scattering measure- 
ment^ has claimed that the value of 7 for x = 0.62 
is reproducible by a phenomenological theory— includ- 
ing the over-damped magnetic excitations for the ferro- 
magnetic fluctuation. However, the additional Ca sub- 
stitution does not lead to ferromagnetism but evolves 
the short-range antiferromagnctic correlation in 0.2 < 
x < 0.5 and the long-range antifcrromagnetic correlation 
in < x < 0.2*£ This experimental fact suggests that 
Ca2- 2; Sr 2 ;Ru04 around x = 0.5 can not be regarded as 
a simple nearly ferromagnetic metal. Furthermore, the 
value of 7o for x = 2 is insensitive to the magnetic fields 
up to 14 T, indicating that the mass enhancement for 
x = 2 is not due to spin fluctuations^ It is thus needed 
to study the role of the ferromagnetic spin fluctuation. 

V. SUMMARY 

In order to clarify the origin of the HF behavior 
around x = 0.5, we have studied the electronic states 
for Ca2- 2: Sr 2 ,Ru04 in 0.5 < x < 2 within the GA on the 
basis of the three-orbital Hubbard model for the Ru t2 g 
orbitals. We have assumed that the Ca substitution af- 



fects the electronic structures mainly by the changes of 
the dp hybridizations between the Ru 4d and O 2p or- 
bitals, and have estimated the mass enhancement on the 
basis of the models taking account of these effects on the 
electronic structures. In particular, we have numerically 
calculated the RF within the GA for three cases with the 
effective models of x — 2 and 0.5 and the special model. 
We have found that the inverse of the total RF becomes 
the largest for the case of x — 0.5, and that the vHs 
for the d xy orbital plays a secondary role in enhancing 
the effective mass. Our results can reproduce the ex- 
perimentally observed tendency of the effective mass in 
0.5 < x < 2 pi 3 - the coefficient of the electronic specific 
heat monotonically increases towards x — 0.5 although 
the vHs appears on the Fermi level at smaller Sr concen- 
tration than x = 0.5. 

Our calculation suggests that the HF behavior around 
x = 0.5 comes from the cooperative effects between mod- 
erately strong Coulomb interaction compared to the total 
bandwidth and the modification of the electronic struc- 
tures due to the rotation of RuOg octahedra (i.e., the 
variation of the dpir hybridizations and the downward 
shift for the d xy orbital). We propose that moderately 
strong electron correlation and the orbital-dependent 
modifications of the electronic structures due to the lat- 
tice distortions play important roles in the electronic 
states for Ca2- :E Sr. r Ru04. 
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